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Abstract We present the results of numerical studies for selected materials with 
strongly correlated electrons using a combination of the local-density approximation 
and dynamical mean-field theory (DMFT). For the solution of the DMFT equations 
a continuous-time quantum Monte-Carlo algorithm was employed. All simulations 
were performed on the supercomputer HLRB II at the Leibniz Rechenzentrum in 
Munich. Specifically we have analyzed the pressure induced metal-insulator tran- 
sitions in Fe203 and NiS2, the charge susceptibility of the fluctuating-valence el- 
emental metal Yb, and the spectral properties of a covalent band-insulator model 
which includes local electronic correlations. 



1 Introduction 

The basic concepts of solid-state physics explain the physical properties of numer- 
ous materials such as simple metals, semiconductors, and insulators. For materials 
with open d and / shells, where electrons occupy narrow orbitals, the additional 
understanding of the role of electron-electron interactions is crucial. In transition 
metals such as vanadium or iron, and also in their their oxides, electrons experience 
a strong Coulomb repulsion because of the spatial confinement in their respective 
orbitals. Such strongly interacting or "correlated" electrons cannot be described as 
embedded in a static mean field generated by the other electrons [1, 2, 3]. The d 
and / electrons have internal degrees of freedom (spin, charge, orbital moment) 
whose interplay leads to many remarkable ordering phenomena at low tempera- 
tures. As a consequence of the competition between different ordering phenomena, 
strongly correlated electron systems are very sensitive to small changes in their con- 
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trol parameters (temperature, pressure, doping, etc.), resulting in strongly nonlinear 
responses, and tendencies to phase separate or to form complex patterns in chemi- 
cally inhomogeneous situations. 

Understanding the metal-insulator transition (MIT) in strongly correlated elec- 
tron systems has been one of the key topics in condensed matter theory since the 
1930's. It was realized early on that the local Coulomb repulsion plays a decisive 
role in the physics of an important class of materials known as Mott insulators, ex- 
amples of which are many transition-metal oxides. While the role of correlations 
due to the on-site Coulomb repulsion has been appreciated for over half a century 
only quite recently the theory progressed to the point where they can be treated 
quantitatively. 

In the last two decades, a new approach to electronic lattice models, the dynami- 
cal mean-field theory (DMFT), has led to new analytical and numerical techniques 
to study correlated electronic systems [1, 3]. This theory - initiated by Metzner 
and VoIIhardt in 1989 - is exact in the limit of infinite dimensions {d ~ °°) [4]. In 
this limit, the lattice problem reduces to a single-impurity Anderson model with 
a self-consistency condition [5, 6]. After the initial studies of conceptually simple 
models, DMFT has advanced to models for real materials and only quite recently 
to purely numerical Hamiltonians obtained from bandstructure calculations, an ap- 
proach which goes by the name LDAh-DMFT. 



2 Computational Method 

2.1 The LDA+DMFT Approach 

In the LDAh-DMFT approach [7] the LDA band structure, represented by a one- 
particle Hamiltonian H^^)^, is supplemented with the local Coulomb repulsion, 
while the correction for double-counting the effect of the local interaction is ab- 
sorbed in /^LDA' 

^ = ^LDA+ H (1) 
m,m' i 

Here, / denotes a lattice site; m and m' enumerate different orbitals on the same 
lattice site. In the present implementation we use an approximate form of the local 
interaction consisting of products of the occupation operators «„„. 

During the last ten years, DMFT has proved to be a successful approach for 
investigating Hamiltonians with local interactions as in Eq. (1)[3]. DMFT treats 
the local dynamics exactly while neglecting the non-local correlations. In this non- 
perturbative approach the lattice problem is mapped onto an effective single-site 
problem which has to be solved self-consistently together with the momentum- 
integrated Dyson equation connecting the self energy E and the Green function 
G at frequency (o: 
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Gr„„A(0)=^J d'k[ico + ^}l-Hl^^{k)-E{(o)]J^,. (2) 

Here, 1 is the unit matrix, jj. the chemical potential, and E{(o) denotes the self- 
energy matrix which is nonzero only between the interacting orbitals. [. . .] ^ ' implies 
a matrix inversion in the space with orbital indices m,m'. The integration extends 
over the Brillouin zone with volume Vb- 



2.2 QMC Method 

The single-site impurity problem can be fomulated without explicit construction 
of the fermionic bath using the effective action formalism [5, 6]. The effect of the 
bath is represented by the frequency-dependent hybridization function A{(o) given 
implicitly by 

G;„\(o) = (o + ii- £,„ - A,„{(o) - I.,„{(o), (3) 
leading to the action 

+ L ^mm'K(T)rm(T)vC'(T)rm'(T) 
+ f dX f dx'Y^rjT)^r„{T-T')W„X^'). (4) 

Jo Jo " 

and the local Green function is obtained by a functional integral over the Grassmann 
variables and i//* 

G,„(T- t') = nwmrWnMwTni^'y^''-''*-^^- (5) 

Here, ^ = j ^^[v^]^[v/'*]exp(5[v/', is the partition function. Such a fomu- 

lation in terms of an effective action in imaginary time is well suited for quantum 
Monte-Carlo (QMC) methods. In the present work we have used predominantly 
the continuous-time hybridization expansion QMC algorithm [8], which consists 
in expanding the exponential in (5) in powers of the hybridization function A and 
sampling the contributions using a QMC random walk, on which almost all of the 
required CPU time is spent upon during the iterative solution of the DMFT equa- 
tions. 

The computational effort of the QMC algorithm based on the expansion of the 
effective impurity action in the impurity-bath hybridization scales as the cube of 
the matrix size [9], which is determined by the mean value of the order of the 
hybridization expansion. The required order of the expansion scales linearly with 
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inverse temperature, and even decreases upon an increase in the electron-electron 
interaction strength and is thus an excellent method for the regime of strong corre- 
lations. Concerning the cluster extension of DMFT, however, the effort of the full 
matrix code grows exponentially with the cluster size. Therefore the hybridization 
expansion solver is efficient for small clusters, but calculations for clusters larger 
than four sites are prohibitively costly. 

From a numerical perspective the continuous-time QMC is perfectly suited for 
runs on the HLRB II. It is almost parallel by definition, and delayed update op- 
timizations [10] are irrelevant to the code since rank-one updates, i.e. single spin 
flips, are not required. Most of the CPU time is spent on the calculation of a trace, 
for which the multiplication of small matrices is needed. These matrices fit in the 
cache and thus the code scales basically linearly with the number of CPUs used. 

The linear scaling in theory is slightly affected in practice by a small overhead 
due to the thermalization sweeps, which need to be performed on each CPU sep- 
arately, and more importantly by the data collection procedure of the ALPS tools 
[11], if more than around 300 CPUs are involved. Our typical jobs, however, are 
rather small scale jobs with 32 to 256 CPUs and CPU times of 10 to 100 CPU hours 
per DMFT iteration, with a typical number of 20 DMFT iterations. The performance 
measured for multiorbital single-site calculations is around 100 to 200 MFlops/core, 
while it fluctuates between 100 and 2000 MFlops/core in simulations of small clus- 
ters using the full matrix code. In a typical scaling test for a 2-site cluster DMFT the 
total CPU time varies between 27 CPU-h on 32 CPUs and 31 CPU-h on 128 CPUs. 
At the same time the performance slightly decreases from 128 MFlops/core on 32 
CPUs to 120 MFlops/core on 128 CPUs. Thus for a typical job, the run time and the 
performance per core scale almost linearly with the number of CPUs, with a slight 
decrease of around ten percent upon increasing the number of CPUs by a factor of 
four. 

The computational requirements of the continuous-time QMC are rather simple. 
It basically needs fast processors in parallel. Not a lot of memory is needed for the 
QMC and, since the random walks on different CPUs are independent of each other, 
fast network communications are not required either. 



3 Results and Discussion 

3.1 Pressure-Driven Metal-Insulator Transition in Hematite 

An important example for of the MIT is the pressure driven transition seen in MnO 
[15], BiFe03 [16] or Fe203 [14], which is accompanied by a change of the local 
spin state (high spin (HS) to low spin (LS) transition). Understanding the pressure- 
driven HS-LS transition and its relationship to the MIT and structural and/or volume 
changes is in fact relevant to a broader class of oxides, often with geophysical im- 
plications. 



Material-Specific Investigations of CoiTelated Electron Systems 



5 



We have studied the spin transition and the MIT in hematite (a-Fe203) under 
pressure using the LDA+DMFT approach [7] including the effects of tempera- 
ture and magnetic long-range order (LRO). At ambient conditions, hematite has 
the corundum structure and is an antiferromagnetic (AFM) insulator with a Neel 
temperature Tn=956 K [17]. The iron ions have a formal Fe^+ valence with five d 
electrons giving rise to a local HS state. Photoemission spectroscopy (PES) clas- 
sified hematite as a charge-transfer insulator [18, 12]. A charge gap of 2.0-2.7 eV 
was inferred from dc conductivity data [19]. Under pressure, a first-order phase 
transition is observed at approximately 50 GPa at which the specific volume de- 
creases by almost 10% and the crystal symmetry is reduced (to the RhoOs-II struc- 
ture) [20, 21, 22]. The high-pressure phase is characterized by metallic conductivity 
and the absence of both magnetic LRO and the HS local moment [20]. Badro et 
al. showed that the structural transition actually precedes the electronic transition, 
which is, nevertheless, accompanied by a sizable reduction of the bond lengths [23]. 

The calculations started with the construction of Wannier based Hamiltonians 
from non-magnetic LDA bandstructures for various specfic volumes spanning the 




Fig. 1 Pressure evolution of the paramagnetic single-particle spectra (T=580 K). The high-spin 
solutions are shown in the left and the low-spin solutions in the right panel. All curves are nor- 
malized to unity, and v denotes the volume relative to the volume at ambient pressure. Left inset: 
comparison of experimental PES [12] and IPS [13] data (points) to the li-spectra (solid line) with 
0.6 eV Gaussian broadening which is comparable to the experimental resolution. Right inset: com- 
parison of the v=0.8 DMFT (smooth curves) and LDA spectra (jagged curves): total rf-electron 
spectra (black), the total oxygen p-electron spectra (blue). Figure taken from Ref. [14]. 
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range of experimentally accessible pressures. Given the nominal configuration 
of the Fe-d shell in hematite, a transition between the high-spin 5 = 5/2 state and 
the low-spin S =1/2 state is expected. Our LDAh-DMFT calculations were first per- 
fomed in the paramagnetic phase and exhibited indeed such a HS-LS transition ob- 
served as a discontinuous drop of the expectation value (5^). The HS-LS transition 
is accompanied by the disappearance of the charge gap and a substantial change in 
the single-particle spectrum (Fig. 1). The comparison to iso-electronic MnO, which 
exhibits a similar HS-LS/insulator-metal transition, reveals a striking difference be- 
tween the two materials. While the transition in MnO sets on before the charge gap 
is closed due to the pressure-induced increase in bandwidth and crystal-field split- 
ting, the transition in hematite is instead characterized by a continuous closing of 
the gap followed by an abrupt spin transition. Comparing these two materials we 
have identified the microscopic mechanism behind the observed transitions [14]. 

A further issue addressed the question of long-range magnetic order. The oc- 
curence of magnetic order in strongly correlated oxides is important also in the 
broader context of first-principles electronic structure methods such as LDA. In sev- 
eral such materials a gap appears in the anti-ferromagnetic LDA solutions, while the 
non-magnetic solutions are metallic. This is sometimes interpreted as antiferromag- 
netism causing the opening of the gap. Our results explicitly showed that magnetic 




Energy (eV) 



Fig. 2 Spin-polarized Fe-d spectra at ambient pressure for 580 K, 1160 K, and 1450 K (top to 
bottom). Left inset: the same spectra averaged over spin. Right inset: staggered magnetization vs. 
temperature for v=1.0 (circles; black), v=0.9 (squares; red), v=0.85 (diamonds; green), and v=0.8 
(triangles; blue), where v denotes the volume relative to the volume at ambient pressure. For empty 
symbols only the low-spin solution exists. The lines represent mean-field fits. Figure taken from 
Ref. [14]. 
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Fig. 3 (a) Charge fluctuations quantified by the mean-square deviation of the / occuancy in Yb, 
Ce, Pr, and Nd at 630 K. If not shown error bars are smaller than the symbol size, (b) Equal-time 
charge fluctuations (diamonds, the same as in left panel) compared to their imaginary-time average 
given by the product of temperature and the static susceptibility Tx^^^ (squares). Figure taken from 
Ref. [24]. 

order has in fact only a marginal effect on the single-particle spectrum and the size 
of the charge gap does not change between the paramagnetic and the antiferromag- 
netic phase (see Fig. 2). We also clarified that the presence of antiferromagentic 
order has only a marginal effect on the HS-LS/insulator-metal transition [14]. 



3.2 Fluctuating Valence and Valence Transition of Yb under 
Pressure 

The valence state of rare-earth atoms in lanthanide compounds also affects their 
physical properties. The determination of the lanthanide valence from first princi- 
ples and the description of the 4/ electrons have been long-standing challenges due 
to the duality between their atomic character, with a strong local electron-electron 
interaction, and their itinerant character due to the lattice periodicity. Theories based 
on a two-species picture, which treat part of the 4/ electrons as atomic and the rest 
as itinerant, succeeded in reproducing the trends across the lanthanide series for 
compounds with integer valence [25] . Nevertheless, besides being conceptually un- 
satisfactory, the two-species picture cannot describe transitions between different 
valence states as well as the heavy-fermion behavior of the charge carriers. 
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We have studied the valence transition in the elemental Yb metal under applied 
pressure. Yb and Eu in their elemental form behave quite differently from the other 
lanthanides. If we define the valence as the number of electrons participating in 
bonding, the majority of the lanthanide series is trivalent, however for Yb and Eu 
the 3+ and 2+ valence states are close to degeneracy with the 2+ state being more 
stable at ambient conditions [25]. This results in a number of anomalous properties, 
such as a larger molar volume compared to the general trend in the lanthanide series, 
and a lower bulk modulus [26]. The thermal-expansion coefficient of Yb is three 
times larger than for most other lanthanides [27]. 

In accord with the experimental data our calculations reveal a continuous de- 
crease of the /-shell occupancy, which corresponds to a crossover from the /'^ to 
the /'^ local charge state as was verified by calculating the single-particle spec- 
tra [24]. The change of the / occupancy is connected to a transfer of electrons 
to spd-hands. The interesting question arises what the observable differences are 
between such self-doping in Yb and the rare-earth materials without valence-state 
degeneracy. The answer is shown in Fig. 3 where the mean-square deviations of the 
/-occupancy and the static local charge susceptibilites as a function of specific vol- 
ume are compared for Ce, Nd, Pr, and Yb. The former three exhibit a monotonous 
increase of the charge fluctuations with pressure and a small charge susceptibility. 
The origin of the increasing charge fluctuations is the growing frequency of hop- 
ping processes involving the / electrons. However, these processes are energetically 
costly due to the on-site Coulomb repulsion and thus very short lived, which leads 
to the small charge susceptibility. The behavior of Yb is strikingly different. Besides 
the increase of charge fluctuations at the smallest volumes, which has essentially the 
same origin as in the other materials, the charge susceptibility vs. volume exhibits 
a local maximum approximately in the middle of the valence-transition region [24]. 
This feature is directly related to the quasi-degeneracy of the /'^ and the /'^ states. 
As a consequence, the charge fluctuations are energetically favorable and thus lead 
to a large susceptibility. The observation of the large charge susceptibility provides 
an explanation for the reported softness of Yb in the valence transition region. 

We conclude that the observed behavior is a common and distiguishing feature of 
fluctuating valence systems characterized by the quasi-degeneracy of local charge 
states. 



3.3 Metal-Insulator Transition in NiS2-xSex 

The NiS2-.rSe;i series presents an important model system in which a metal-insulator 
transition can be controlled either by varying the Se content x, temperature T, or 
pressure P [30, 31, 32, 33, 34]. Despite a vast amount of available experimental 
data a satisfactory material-specific theory for the microscopic origin of the MIT in 
NiS2-j:Sejc is still missing. Using the LDAh-DMFT approach we have discovered an 
unexpected mechanism which controls the opening of the charge gap in NiS2-;tSeA^. 
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NiX2 (^=S, Se) can be viewed as NiO with the O atom replaced by an X2 dimer 
(Fig. 4), which accommodates two holes in its p*„ anti-bonding orbitals leading to 
an ^2 ^ valence state. It was previously asserted that the empty p'^ bands do not 




Fig. 5 The calculated (U=5 eV, 7=580 K) orbitally resolved spectra of NiS2 [Ni-f/:e^ (black) 
and t2g=e^ (thick red) + a\g (thin red); S-p: total (blue) and (blue-shaded)] compared to the 
experimental XPS-fBIS data (gray, shaded) [28]. The inset shows the conduction band spectra 
for U=5, 6, 7, and 8 eV. Figure taken from Ref. [29]. 
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play an active role in the physics of NiX'2. If true the analogy to NiO becomes com- 
plete in the sense that a charge-transfer gap forms between the ligand /9-band and 
the upper Ni-^f Hubbard band, whose position is sensitive to the strength of the on- 
site Coulomb repulsion U . Within such a scenario a stronger screening leading to 
smaller U in NiSe2 can explain why NiS2 is an insulator while NiSe2 is a metal. 
However, our calculations revealed only a moderate difference in the Coulomb re- 
pulsion between NiS2 and NiSe2, and their respective groundstates turned out to be 
unaffected by variations of U . This is readily understood by the observation that it is 
not the upper Ni-c/ Hubbard band, but rather the X-p% band, which forms the bottom 
of the conduction band (see Fig. 5). In the meantime this finding was corroborated 
by X-ray emission and absorption measurements at the K-edge of sulfur [29]. 

Our calculations provide the following picture. The different groundstates of 
NiS2 and NiSe2 are consequence of a larger bonding-anti-bonding splitting within 
the S-S dimer The splitting in the Se-Se dimer is too small to open a gap. The 
closing of the gap in NiS2 by the application of pressure is due to a broadening of 
the bands while the S-S bond-length related bonding-anti-bonding splitting remains 
unchanged since the S-S dimer behaves as a rigid object in a soft matrix. 



3.4 Interaction Driven Insulator-to-Insulator Transition 

While the Hubbard model has become a paradigm for the description of electronic 
correlations in metals and the metal-insulator transition [2], much less attention 
has so far been paid to electronic correlations in band insulators. E.g., modeling 
of Kondo insulators is far from trivial, and the only recently achieved progress in 
the topological classification of band insulators [35] demonstrates that our under- 
standing of the insulating state is still incomplete. 

Motivated by several investigations of the ionic Hubbard model [37, 38, 39, 40, 
41, 42, 43, 44] we have analyzed a covalent insulator as a complementary example 
of a band insulator As a covalent insulator we denote a band insulator with partially 
filled local orbitals. This definition implies that the band gap is a hybridization gap 
arising from a particular pattern of hopping integrals. It has been proposed that sim- 
ilar characteristics apply to materials such as FeSi, FeSb2 or CoTiSb [45], some of 
which exhibit temperature dependent magnetic and transport properties reminiscent 
of Kondo insulators. 

In our model study we use a simple particle-hole symmetric model at half-filling 
described by the Hamiltonian 



H 



ka \ ka / ia 



ni]ani[ai 



(6) 




(7) 
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Fig. 6 The T-U phase diagram at fixed V = 0.5 shows the band and Mott insulating phases with a 
first-order phase transition upon increasing the interaction strength U below a critical temperature. 
In the coexistence region both band and Mott insulating solutions of the DMIT equations are found 
depending on the initial guess for the self-energy. For temperatures above the critical end point of 
the coexistence region, there is a regime where the spectral function has a single peak at the Fermi 
energy accompanied by broad Hubbard bands. All energies are given in units of f = 1. Figure taken 
fiom [36]. 



with two semi-circular electronic bands of widths 4r (r = 1 in the following) and dis- 
persions Ek and — e^,, respectively, corresponding to two sublattices coupled by the 
ik-independent hybridization V and a local electron-electron interaction of strength 
U . Here = ^Ja^ia measures the number of electrons with spin a =|, J, on site 
i of sublattice a = a,b. 

We use DMFT in conjunction with the recently developed continuous-time QMC 
[8] to study the evolution from the band insulator at small to the Mott insulator at 
large Coulomb interaction strength. The insulator-insulator transition is discontin- 
uous at finite but low temperatures, with a region of interaction strengths where 
hysteretic behavior with two solutions of the DMFT equations is observed (see Fig. 
6). The behavior of charge and spin gaps upon increasing the interaction strength 
is shown in Fig. 7 [36]. Surprisingly we find that both gaps shrink with increas- 
ing Coulomb repulsion. This behavior is in contrast to the correlation-induced Mott 
insulator where the charge gap increases with increasing interaction strength. Fur- 
thermore, in the correlated insulator charge and spin gaps deviate from each other, 
and the spin gap is smaller than the charge gap. From the self-energy we extract 
a renormalization factor Z defined in the same way as the quasiparticle weight in 
a Fermi liquid. In the band insulator the Z factor describes the renormalization of 
the charge gap at moderate interaction strengths (see Fig. 7). Thus we obtain the re- 
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Fig. 7 Spin and charge gaps in the correlated hand insulator as a function of U for y = 0.5 as 
determined from the spectral function and the spin susceptibility at T = 1 /30, respectively. The 
thick solid line shows the charge gap obtained from second order perturbation theory. The squares 
represent the Z factor Inset: Discontinuous change of spin and charge gaps at the band insulator to 
Mott insulator transition. In the Mott insulator the spin gap 4., = 0. All energies are given in units 
of / = 1. Figure taken from [36]. 

markable finding tliat a concept from Fermi liquid theory can be applied to quantify 
correlation effects in a band insulator! 



4 Conclusions 

We have presented several examples for the application of dynamical mean-field 
theory to investigate quantitatively the properties of selected electronically corre- 
lated materials. As the main numerical tool a continuous-time quantum Monte- 
Carlo algorithm was used to solve the auxiliary impurity problem in the DMFT 
self-consistency cycle. The present method proved to be a powerful tool to calculate 
spectral properties of materials with strong electronic correlations, as well as phase 
transitions near and above room temperature and at arbitrary pressure. The presented 
model analysis for a covalent band insulator with local Hubbard-type interaction 
provided new insights into the previously unexplored effects of correlations in band 
insulators. 
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